Lepton interferometry in relativistic heavy ion collisions - a case study 
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We propose intensity interferometry with identical lepton pairs as an efficient tool for the estima- 
tion of the source size of the expanding hot zone produced in relativistic heavy ion collisions. This 
can act as a complementary method to two photon interferometry. The correlation function of two 
electrons with the same helicity has been evaluated for RHIC energies. The thermal shift of the p 
meson mass has negligible effects on the HBT radii. 
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I. INTRODUCTION 

In high energy heavy ion collisions two-particle 
correlations have been extensively studied both ex- 
perimentally [1 H Q and theoretically 0, II @ 
to obtain direct information about the size, shape 
and dynamics of the source at freeze-out. This 
is usually done via selection of transverse momen- 
tum and rapidity of the correlated particles. Such 
calculations are based on the fact that identical 
particles occurring nearby in phase-space experi- 
ence quantum statistical effects resulting from the 
(anti)symmetrization of the multi-particle fermion 
(boson) wave function. For bosons, the two-particle 
coincidence rate shows an enhancement at small mo- 
mentum difference between the particles and an op- 
posite behaviour is observed for fermions. The mo- 
mentum range of this enhancement or depletion can 
be related to the size of the particle source in co- 
ordinate space. 

Recently, it has been argued that Hanbury-Brown 
and Twiss (HBT) interferometry is a sensitive probe 
of the QCD equation of state and hence formation 
of QGP 0,0 111- Further, it has been argued that 
in contrast to hadron s, two-p article intensity inter- 
ferometry of photons [l(J, fill Ilia [jjj which are pro- 
duced throughout the space-time evolution of the re- 
action and which suffer almost no interactions with 
the surrounding medium can provide information on 
the the history of the evolution of the hot matter 
created in heavy ion collisions |l 1| . 

The two fermion interferometry is a well known 
method used in nuclear physics to estimate the 
size of the source formed after nuclear colli- 
sions ^5|( see [13 f° r a review). In case of fermion in- 
terferometry the symmetric (anti-symmetric) space 
part of the wave function is coupled with the anti- 



symmetric (symmetric) spin part of the wave func- 
tion. In this work we will evaluate correlation func- 
tion for the symmetric spin part and hence anti- 
symmetric space wave functions. Statistical spin 
mixture has been neglected here |l6j . 

In our previous work on photon interferome- 
try the spin-averaged source function was used 
to study the two photon correlation functions. Here 
we use spin-dependent source functions to evalu- 
ate the two-lepton correlation function. The main 
aim of this study is to extract the HBT radii of 
the source from spin-dependent electron correlation 
function and show that the HBT radii extracted here 
is similar to those obtained from the two photon in- 
terferometry. This indicates that the results from 
photon interferometry is not very sensitive to the 
spin of the photon. In addition to this it will be 
useful to see the effects of the shift of vector meson 
masses with temperature on the correlation function 
of two electrons. 

The paper is organized as follows. In the next sec- 
tion we give a general discussion on the correlation 
function for fermions and the associated kinematics. 
This is followed by the section which deals with the 
space time evolution. Section IV, deals with the re- 
sults on two lepton correlations at RHIC energies, 
results with mass variation of vector mesons in a 
hot medium and brief discussion on the results of 
two electron interferometry vis-a-vis two photon in- 
terferometry. We summarize our findings in the last 
section. 



II. CORRELATION FUNCTION 

The two-particle correlation function in momen- 
tum space is defined as, 
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where k\ and k 2 are the three momenta of the two 
particles. P\(ki) and P 2 (ki,k 2 ) represent the one- 
and two- particle inclusive electron spectra respec- 
tively. These are defined as, 



and 



P x (k) = / d 4 x S(x,k) 



(2) 



P2(ki,k 2 )=Pi(k 1 )P 1 {k 2 ) 



d i x 1 d 4 x 2 S(xx,K) S{x 2 ,K) cos(Aa; M Afc / ,) 



(3) 



where, K = {k\ + fc 2 )/2, Afc^ = k lf2 - k 2ll = q^, x { 
and ki are the four-coordinates of position and mo- 
mentum respectively and S(x, k) is the source func- 
tion, which defines the average number of electrons 
with four-momentum k emitted from a source el- 
ement centered at the space-time point x. In the 
present case S(x,k) is the thermal emission rate of 
electrons per unit four volume. For processes of the 
form a(pi)+ /3(p 2 ) — > e + (k\) + e~(fc), we have |l7j . 



S(x, k) 



ds 
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, (4) 
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where a and (3 are either quarks or pions, ki = 
(Ei,ki) is the four vector for the particle i. The 
masses of quarks and electrons are neglected here. 
TV is the overall degeneracy of the particles a and 
(3, f(Ei) denotes the thermal distribution functions 
and s, t, u are the usual Mandelstam variables. The 
expressions for a, b, c, and the integration limits, 
Eimin, E 2mm and E 2max are given in Ref. [T3| - 
The spin-dependent invariant amplitude, Ad, for 
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the processes q R q~l — > 

e R e l( e L e R) and 
been calculated using standard field theoretic tech- 
niques [lgj . 

We shall be presenting the results as a function of 
longitudinal (qiong), outward (q ou t), sideward (q s ide) 
and invariant momentum differences {qi nv ) of the 
two leptons and these are defined as, 



qiong — k\ z — k 2z 

= kir sinh(yi) - k 2T sinh(y 2 ) 



(5) 
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= (7) 



= -2fc 1T fc 
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cosh(yi - y 2 ) - cos(-0i - ^2) 
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where, <7r = k\x — k 2 T, Kt — (kiT + k 2 T)/2 with the 
subscript T indicating the transverse component, yi 
is the rapidity and tpiS are the angles made by kiT 
with the x-axis. 

The HBT radii can be extracted by parametrizing 
the calculated correlation function in the following 
Gaussian form: 
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(9) 

The superscript l Asym' indicates that antisymmet- 
ric space wave functions have been used. Also note 
that a negative sign appears before the exponential 
(as opposed to a positive sign in case of the two- 
boson correlation function). The A in Eq[5]is called 
chaoticity parameter, A = for a coherence source 
and A = 1 for a completely chaotic source. We have 
set A = 1 as we are concerned here with a completely 
thermalized system from where the leptons are emit- 
ted. If the source function has more than one max- 
ima then Gaussian parametrization fails to describe 
the structure of the correlation function at large val- 
ues of outward, sideward or invariant momenta de- 
fined in section II. For "binary" source functions like 
mixed phase in the present context, more than one 
maxima in correlation functions are possible G] . The 
values of C^^" 1 will lie between and 1 for a chaotic 
source. 

A gross idea of the source size can also be obtained 
from Ri nv which is defined as: 
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III. SPACE-TIME EVOLUTION 

For the evaluation of two-electron correlations at 
RHIC energies we will consider a scenario where 
QGP is formed in the initial state which evolves 
with time into a hadronic phase via an intermedi- 
ate mixed phase in a first order phase transition sce- 
nario. The mixed phase is a mixture of both quark 
matter and hadronic matter, with the fraction of 
quark matter decreasing with time to zero when the 
phase transition is complete. The hot hadronic gas 
then expands till the system freezes-out. We also 
evaluate the correlation function when the vector 
meson masses {p in this case) varies with tempera- 
ture according to the universal scaling scenario pro- 
posed by Brown and Rho (BR) |l9|. 

The initial condition for RHIC energies in terms 
of the initial temperature (Tj) is calculated from the 
number of particles per unit rapidity at the mid- 
rapidity region for those energies according to the 
following equation, 

3 2tt 4 1 dN _ 1 dS ( \ 

where dS(dN) is the entropy (number) contained 
within a volume element AV = it R 2 ti dy, R as 
the radius of the colliding nuclei. g e ff is the effec- 
tive statistical degeneracy, £(3) denotes the Riemann 
zeta function and y is the rapidity. For massless 
bosons (fermions) the ratio of dS to dN is given by 
2?r 4 /(45C(3)) ~ 3.6 (4.2), which is a crude approx- 
imation for heavy particles. For example the above 
ratio is 3.6 (7.5) for 140 MeV pions (938 MeV pro- 
tons) at a temperature of 200 MeV. — 7r 2 t/fc/90 
is determined by the statistical degeneracy (gk) of 
the system formed after the collision. Taking the 
particle multiplicity per unit rapidity to be 1100 for 
Au+Au collisions at RHIC, weget T; = 264 MeV 
for an initial time of 0.6 fm/c [2(|. The critical tem- 
perature (T c ) is taken to be 170 MeV here. The 
freeze-out temperature (Tf) is taken to be 120 MeV, 
which describes the transverse momentum distribu- 
tions of produced hadrons |22| with the same space 
time evolution model used here. For the equation of 
state (EOS) which plays a central role in the space 
time evolution, we have considered a hadronic gas 
with particles of mass up to 2.5 GeV and includ- 
ing the effects of non-zero widths of various mesonic 
and hadronic degrees of freedom [2j| . The velocity of 
sound (c s ) corresponding to this EOS at freeze-out 
is about 0.18 [23j. The bag model EOS has been 
used for the QGP phase. Using the above inputs 
and assuming boost invariance along the longitudi- 
nal direction [24( the (3+l)-dimensional hydrody- 
namic equations have been solved to study the space 
time evolution [2{| from the initial QGP phase to 
freeze-out with an intermediate mixed phase of QGP 
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FIG. 1: Transverse momentum distribution of single 
electron for Au + Au collisions at 200 AGeV at RHIC. 
Dotted (dashed) line indicate results for QGP (Mixed) 
phase and dot-dashed (solid) line represent pr spectra 
for hadronic (total) phase. 



and hadrons. The initial energy density profile used 
here is same as the one used in |2Gj and the initial 
radial velocity is taken as zero. 



IV. RESULTS 

In this section we present the results of two- 
electron interferometry at RHIC energies. In Fig. ^ 
the transverse momentum distribution of the elec- 
trons for Au -I- Au collisions are depicted for various 
phases. The initial temperature is taken as Tj = 264 
MeV. We have assumed that the mass of the inter- 
mediary p in pion annihilation varies according to 
Brwon-Rho scaling [19|. 

The results on the interferometry can be presented 
for several combinations in the variables y, tp and 
hr- For simplicity, we will present all two-lepton 
correlation functions for single electrons with mo- 
mentum 1 GeV/c, with the assumption that around 
this value of transverse momentum most of the elec- 
trons will have a thermal origin. We take ^2 = 
and 2/2 = for all cases varying ipi and y\ wher- 
ever necessary. Figs. EJa-d) show the variation of 
the correlation strength (C^ sym ) as a function of 
qiong, (lout, Iside and q lnv for various phases (sum = 
QGP+Mixed+hadrons). The HBT dimensions ex- 
tracted from these correlation functions are shown in 
Table I. To obtain the longitudinal dimension Ri on g 
in the longitudinally co-moving system (LCMS) of 
reference one should multiply the numbers given in 
Table I by the Lorentz factor (= cosliyx), where 
yx is the rapidity corresponding to K defined in 
section II. The width of the correlation function for 
the hadronic phase is the largest as compared to the 
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FIG. 2: Correlation function, C , ^ sym as a function of 
qio ng , qout, qside and q inv for Au + Au collisions at 200 
AGeV at RHIC. Solid (dashed) line indicate results for 
QGP (Mixed) phase and dot-dashed (dotted) line repre- 
sent correlation function for hadronic (sum) phase. 
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FIG. 3: Same as Fig. EJb-c) when the mass of the p 
varies with temperature according to BR scaling. 



other phases along q ou t, followed by that for the 
mixed phase and the QGP phase. For qw e , the 
values are comparable for QGP and mixed phases, 
while it is slightly lower for hadronic phase. The 
HBT dimensions satisfy the relation R ou t /Rside ~ 2 
for the correlation functions denoted by 'sum' in the 
Table I. This is consistent with the earlier calcu- 
lations on pion interferometry j^l that the ratio 
Rout I Rside will be larger than unity in a first order 
phase transition scenario due to the appearance of 
mixed phase and hence time delay due to the slow 
down of the expansion rate. We would like to empha- 
size here that the HBT radii give the length of ho- 
mogeneity of the source (2^| and this is equal to the 
geometric size if the source is static. However, for a 
dynamic source e.g. the system formed after ultra- 
relativistic heavy ion collisions, the HBT radii is 
smaller than the geometric size (see 29, 30, 31, 32]). 

As indicated earlier, we also consider a scenario 
where the p meson mass varies with the tempera- 
ture of the medium. The obvious motivation is to 
comment on this very important issue based on two- 
electron interferometry. In Figs. |3fa-b) we show the 
variation of the correlation strength (C2) as a func- 
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tion of q ou t and q S ide for various phases with the 
masses of the p modified in the medium according 
to BR scaling. Unfortunately, these results are not 
very different from the corresponding ones without 
medium effects shown in Fig. We do not present 
the results for qi ong and qi nv as they are very similar 
to the case without mass variation. All these indi- 
cate that lepton interferometry is not, probably, a 
suitable probe to detect the in-medium modification 
of vector mesons. 

We also note that the HBT radii of various 
phases extracted from photon interferometry |10| 
with spin-averaged source function and those ob- 
tained here from electron interferometry with the 
spin-dependent source functions are similar. This 
indicates that the results from photon interferome- 
try is not very sensitive to the spin of the photon. 
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Table 1 : Values of the various HBT radii in fm. 



V. SUMMARY AND OUTLOOK 

The two-electron correlation functions have been 
evaluated for RHIC energies. The spin dependence 
of the source functions of the electrons originating 
from the QGP phase and hadronic phase have been 
considered explicitly through the invariant ampli- 
tude. (3+1) dimensional relativistic hydrodynamics 
has been used for the evolution in space and time 
from the initial QGP phase to the final hadronic 
phase with intermediate mixed phase of QGP and 
hadrons in a first order phase transition scenario. 
In contrast to the invariant mass distribution of the 
dileptons [HEH, the HBT radii obtained from the 
two electrons interferometry is seen to be insensi- 
tive to the in-medium modifications of the inter- 
mediary vector mesons, apparently because of the 
cancellation of such effects between the numerator 



and the denominator of Eq. ^ The results obtained 
from the electron interferometry here is similar to 
those obtained from the photon interferometry pH | 
(where spin averaged amplitudes were used for pho- 
ton production). The values of HBT radii extracted 
from the two-lepton correlation functions show that 
Rout/Rside ~ 2, consistent with the assumption of 
a first order phase transition scenario [2]j. It has 
been checked that the HBT radii extracted from the 
correlation function evaluated by using symmetric 
space and anti-symmetric spin wave function is sim- 
ilar to the values given in Table 1. Experimentally 
it is difficult to verify the predictions made here, 
however, first attempt to measure two photon cor- 
relation functions in heavy ion collisions has already 
been made by WA98 collaboration [lflj . 

Unlike hadrons which are dominantly emitted 
from the freeze-out surface of the fireball, the lep- 
tons and photons are produced and emitted from 
all the evolution stages of the matter formed after 
the nuclear collisions. Therefore, lepton and pho- 
ton interferometry can give, in principle, 'the length 
of homogeneity' of the system at any stage of the 
evolution. We have not come across any theoretical 
work on lepton interferometry with spin dependent 
invariant amplitude and (3+l)dimensional hydrody- 
namical expansion. However, there are scopes of im- 
provements of the present work. A first order phase 
transition is assumed here in absence of satisfactory 
understanding of the order of the phase transition 
from lattice QCD [U|3|. It will be interesting to 
look into the two particle correlation functions with 
a continuous transitions. Near the phase transition 
the quarks and gluons may not behave as massless 
particles due its interactions with other particles in 
the thermal bath although they are treated here as 
massless. The bag model equation of state has been 
used here for simplicity, work with more realistic 
EOS 21] is necessary. We have assumed that e~ 
from annihilation of thermal quarks and pions dom- 
inates at pt = 1 GeV. To justify this a detail analy- 
sis of the pt distribution of e~ 's from other sources, 
e.g. Dalitz decays, open charm decays, Drell-Yan 
processes etc are required. A detailed calculations 
taking in to account some of these factors will be 
published elsewhere. 
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